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Nonzero ground state entropy, 5*0 7^ 0, is an important subject in statistical mechanics as 
an exception to the third law of thermodynamics (e.g., |jl|). This is equivalent to a ground 
state degeneracy per site W > 1, since 5*0 = ksliaW. The g-state Potts antiferromagnet 
(AF) 0,^ exhibits nonzero ground state entropy (without frustration) for sufficiently large 
g on a given lattice A, or more generally, a given graph G, and serves as a valuable model 
for the study of this phenomenon. The zero-temperature partition function of the above- 
mentioned g-state Potts AF on G satisfies Z{G,q,T = 0)paf = P{G,q), where P{G,q) 
is the chromatic polynomial (in q) expressing the number of ways of coloring the vertices 
of the graph G with q colors such that no two adjacent vertices have the same color [§,^. 
Thus, iy({G}, q) = lim„_^oo P{G, qY^"', where n = v{G) is the number of vertices of G [|],0 
and {G} = lim„^ooG. W{{G},q) has been calculated exactly for the triangular lattice 
P) and various families of graphs 0, 1^- ||15|. The special values for the square (sq) and 
kagome (kg) lattices W{sq,3) |T^ and W{kg,3) (which can be extracted from ||T7||8|) are 



also known. However, aside from the triangular case, the exact calculation of W^(A, q) for 
general q on lattices A of dimensionality d > 2 remains an outstanding challenge. In this 
work we report exact calculations of W on infinite-length, finite-width strips of the square 
and kagome lattices that exhibit the analytic properties expected for the W functions on the 
respective full 2D lattices and, in this sense, constitute the closest exact results that one has 
to these W functions. 

Let us describe these analytic properties. Denote lim„^oo C = {G}- Since P{G,q) is a 
polynomial, one can generalize q from Z+ to M and indeed C W{{G},q) is a real analytic 
function for real q down to a minimum value, qc{{G}) |]7|,[T8|. For a given {G}, we denote the 



continuous locus of non-analyticities of W as B. This locus B forms as the accumulation set 
of the zeros of P(G, q) (chromatic zeros oiG) Bsn ^> 00 |lT9|-p|,[^ and satisfies B{q) = B{q*). 



In cases where B serves as a natural boundary, dividing the q plane into different regions, 
W has different analytic forms in these different regions. The point qc is the maximal point 
where B intersects the real axis, which can occur via B crossing this axis or via a line segment 
of B lying along the axis. The chromatic polynomial P{G, q) has a general decomposition 
as P{G,q) = Co(g) + J2j(^j{l){^jil)y^"' where the aj{q) and Cj^o{q) are independent of n, 
while Co(g) may contain n-dependent terms, such as (—1)"", but does not grow with n like 
{const.)"' with |corast.| > 1, and tj is a G-dependent constant. A term a^(g) is "leading" 
{£) if it dominates the n ^ 00 limit of P{G,q). The locus B occurs where there is an 
abrupt nonanalytic change in W as the leading terms a^ changes; thus the locus B is the 
solution to the equation of degeneracy of magnitudes of leading terms. Hence, W is finite 
and continuous, although nonanalytic, across B. 

From exact calculations of VT on a number of families of graphs we have inferred several 



general results on B: (i) for a graph G with well-defined lattice structure, a sufficient condi- 
tion for B to separate the q plane into different regions is that G contains at least one global 
circuit, deffned as a route following a lattice direction which has the topology of the circle 
S^ and a length ig,c. that goes to infinity as n ^ oo |[10|J2^ . For a (i-dimensional lattice 



graph, the existence of global circuits is equivalent to having periodic boundary conditions 
(BC's) in at least one direction. Further, (ii) the general condition for a family {G} to have 



a locus B that is noncompact (unbounded) in the q plane |jTl| shows that a sufficient (not 
necessary) condition for {G} to have a compact, bounded locus B is that it is a regular 
lattice [llll , p!2| , [l5| . The third and fourth general features are that for graphs that (a) contain 
global circuits, (b) cannot be written in the form G = Kp + H P,|23[|, and (c) have compact 
B, we have observed that B (iii) passes through q = and (iv) crosses the positive real axis, 
thereby always defining a qc- 

From exact calculations of W for a number of infinite-length, finite-width (homogeneous) 
strips of 2D lattices with free boundary conditions in the longitudinal direction (and free or 



periodic BC's in the transverse direction) |1T0| , [13| , it is found that the resultant loci B consist 



of arcs (and possible real line segments) which, although compact, do not separate the q 
plane into different regions, do not pass through g = and, for the arcs, do not necessarily 
intersect the real q axis. These calculations showed that as the strip width Ly increases, the 
complex-conjugate (c.c.) arcs comprising B tend to elongate so that the gaps between them 
decrease, and the left endpoints of the c.c. arcs nearest to g = move toward this point, 
thereby leading to the inference that in the limit Ly ^ oo, these arcs will close to form one 
or more regions, and B will pass through g = and will cross the positive real axis at one or 
more points, thereby defining a qc- In turn, this motivates the conclusion that the properties 
(i)-(iv) hold for W{A,q) and ;B on a lattice A (in the thermodynamic limit, independent of 
the boundary conditions used). The advantage of cyclic strip graphs is that these properties 
are present for each finite Ly rather than only being approached in the limit Ly ^ oo as for 
open strips. 

Our method for obtaining exact W functions that exhibit the analytic properties expected 
for W on a 2D lattice is as follows. We calculate P{Ga, q) on L^ x Ly strips of the lattice 
A with periodic (i.e., cyclic) BC's in the longitudinal {L^) direction, then take L^, — >■ cx) and 
calculate W and the resultant B. By construction, these W functions and the associated 
loci B embody the four general properties given above. For each strip, the exterior of B in 
the q plane, denoted as the region i?i, is the maximal region into which one can analytically 
continue W from the real interval q > qc- The calculation of W for a cyclic strip of a given 
width is considerably more difficult and the result more complicated than that for the open 
strip (free L^ BC) of the same width; the value of the cyclic strips is that the resultant W 



exhibits the analytic features of the full 2D W function. The boundary condition in the 
transverse direction is not important for these results since the width is finite; for simplicity 
we use free transverse BC's. 



We use an extension of the generating function method of Ref. |T^ from open to cyclic 
strip graphs Ga- The generating function r{GA,ci,x) yields the chromatic polynomials for 
finite-length strips of A as the coefficients in its Taylor series expansion in the auxiliary 
variable x about x = 0. Here, T{GA,q,x) = M{G\,q,x)/V{G\,q,x), where J\f and V are 
polynomials in x and q (with no common factors). The degrees of these, as polynomials in 
X, are denoted j^ax = deg^iJ^) and k^ax = deg^iV). The J\f are not needed here (they 
will be given elsewhere) since W and B are determined completely by V, independent of A/" 
10| . For a particular G\, writing V = YYjlTi^ ~ Xjx), W is given in region Ri and |iy| in 



other regions ||2^ hj W = {XmaxY and |W^| = |Amax|*, where Xmax denotes the A in P with 
maximal magnitude in the respective region and t = Lx/n = l/Ly for the square strip and 
1/5 for the kagome strip considered here. 

We first consider cyclic strips of the square lattice. For Ly = 1, B consists of the unit 
circle |g — 1| = 1 so that qc = 2, and W = q — 1 for q G Ri. For Ly = 2, from the known P 
||T9| , we found that B separates the q plane into four regions, qc = 2, and W = (g^ — 3g + 3)^/^ 
for q E Ri [0]. We have calculated the generating function for the Ly = 3 case. This has 
jmax = 8 and k^ax = 10 and is considerably more complicated than the Ly = 3 open strip, 
where jmax = 1 and kmax = 2. For V we find 



V{sq{Ly = 3), g, x) = (1 + 6,g,iix + &,g,i2a;^)(l + &.g,2ia; + &sg,22a;^ + ^,233;' 



3-\ 



X 



(1 + x) [1 + (g - 2fx] [1 - (g - 2)x] [1 - (g - 4)x] [1 - (g - l)x] (1) 
where 

&s,,ii = -(9-2)(g'-3g + 5) (2) 

bsg,i2 = (g - 1) (g' - 6g2 + 13g - 1 1) (3) 

bsq,2i = 2g' - 9g + 12 (4) 

bsq,22 = g' - 10g=^ + 36g2 - 56g + 31 (5) 

bsg,23 = -(g - l)(g' - 9q^ + 29g2 - 40g + 22) (6) 
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FIG. 1. Locus B for W for the cx) x 3 cyclic strip of square lattice. Chromatic zeros for L^ — 20 (i.e., 
n = 60) are also shown. 

The boundary B is shown in Fig. |I]. It divides the q plane into several regions and crosses 
the positive real axis at Qc = 2.33654 and q = 2. Thus, this Ly = 3 cyclic strip is the first 
one sufficiently wide as to yield a value of Qc above g = 2; indeed, the value of qc for this 
strip is only about 20 % below the value for the full 2D lattice, viz., gc = 3 0]. In region Ri 
including the real interval q > qc, 



W{{G 



sq{Ly=3) 



},qeRi 



-1/3 



(g-2)(g2-3g + 5) + 



(g^ - 5g + 7)(g^ - 5g^ + llg^ - 12g + 8) 



1/2- 



1/3 



(7) 



At qc, W = 1.18487. In the region that includes the real interval 2 < q < qc, \W\ = |g — 4|^/^. 
In the region that includes the real interval < q < 2 and in the regions centered at 
roughly q = 2.4 ± 0.9i, \W\ is given by the respective maximal cube roots of the equation 
C^ + bsg^2i^^ + bsq,22(,+bsq,23 = 0. As an algebraic curve, B has several multiple points (defined 
as points where several branches of this curve cross intersect). 

We next consider a cyclic strip of the kagome lattice comprised of m hexagons with each 
pair sharing two triangles as adjacent polygons (as in Fig. 1 in [|10| for the open strip). F 



has jmax = 8 and kmax = 9 as compared with jV 



1, kmax = 2 for the open strip of the 



same width 10 . We calculate 



V{kg{Ly = 2), g, x) = (1 + hkg^uX + bkg^uX^) (1 + bkg,2ix + bkg,22X^){l + bkg,3iX + bkg,32X^) x 



where 



[1 - (g - 2)x][l - (g - 4)x][l - {q - l)(g - 2fx] 



Jkg,U 



-{q - 2)(g^ - 6g^ + Uq^ - 16g + 10) 



(8) 



(9) 



bkg^i2 = {q-lf{q-'2f 



(10) 



bkg,2i = -q^ + 7g' - 19g + 20 



(11) 



bkg,22 = {q-l){q-2f 



6fc,,3i = ll-9g + 2g^ 



(12) 
(13) 



fe%,32 = -(g-i)(g-2)' 



(14) 



Define Afcc,j,± = {l/2)[—bkgji ± (^^gji — 46fcc,j2)"'^^^]- Again, B divides the q plane into several 
regions (Fig. 2). In region Ri, W is determined by \kg,i,+: 



W{{G,g^,^=2)},q) = 2-'/'{q-2y/' 



q^ - 6q^ + Uq^ - 16g + 10 



+ 



q^ - I2q^ + 64g^ - 200g^ + 404g^ - 548g^ + 500g^ - 292g + 92 



1/2 



1/5 



(15) 



As is evident from Fig. 2, the value of qc is within about 10 % of the inferred exact value 



gc = 3 for the 2D kagome lattice |]14|. It is impressive that an infinite strip of width L, 



yields a qc this close to the value for the full 2D lattice. 
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FIG. 2. Locus B for W for oo x 2 cyclic strip of the kagome lattice. Chromatic zeros for m — 10 {n — 50) 
are also shown. 



Another interesting feature of these resuhs is the fact that the chromatic zeros and their 
accumulation set B contain support for Re{q) < 0. This is in contrast with the situation for 
strips with free longitudinal BC's |10[ and provides further support for our earlier conjecture 
that a necessary condition for this Re{q) < feature is that the graph family have global 
circuits. 

We have also computed W and B for the cyclic strip of the triangular strip with Ly = 2. 
We find V = {1 — x)[l — {q — 2)^x][l + (2g — 5)x + {q — 2)^x^]. B separates the q plane into 
three regions and crosses the positive real axis at gc = 3 and at g = 2. The qc value for this 
strip is one unit less than the value gc = 4 for the full 2D lattice. 

Similar calculations can be carried out for infinite-length cyclic strips G\ of greater 
widths. Our method can also be applied to lattices with d > 3. To do this, one would 
use the generating function method to calculate P for tubes with longitudinal PBC's and 
successively larger [d — l)-dimensional cross sections. We believe that this application, as 
well as that to other 2D lattices, is promising. 
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